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ABSTRACT 

In  a  previous  paper  devoted  to  the  numerical  solution  of  the  Stefan 
problem,  the  author  has  proposed  a  numerical  scheme  to  solve  the  heat 
equation  on  a  variable  mesh;  this  scheme  is  a  generalization  of  the 
classical  Crank-Nicolson  scheme  since  it  is  identical  to  the  Crank- 
Nicolson  scheme  in  the  particular  case  of  a  fixed  mesh.  Numerical 
experiments  have  been  performed  in  one  and  two  space-dimensions,  but  no 
mathematical  results  had  been  proved.  In  the  present  paper,  the  sta¬ 
bility  and  convergence  of  the  scheme  are  established  together  with  an 
error  estimate. 
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sluNIKl.-ANOl-:  AND  KXl'I.ANAT  1  ON 


Moving  boundary  problems  occur  in  many  areas  of  physics  and  engineer¬ 
in').  The  best  known  example  is  the  Stefan  problem  which  represents  the 
melting  of  a  solid:  the  moving  boundary  is  the  interface  between  the 
solid  and  the  liquid)  in  each  phase,  the  temperature  satisfies  the  heat 
equation.  An  efficient  way  of  solving  moving  boundary  problems  is  to  use 
a  variable  mesh.  In  a  previous  paper,  such  a  method  has  been  projxised 
for  the  heat  equation  and  numerical  experiments  have  been  performed.  In 
this,  pa  pet  ,  the  stability  and  convergence  of  this  method  are  proved. 
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STABILITY  AND  CONVERGENCE  OF  A  GENERALI  I'll »  RANK-ND  •  I..  N 
SCHEME  ON  A  VARIABLE  MESH  FOR  THE  HEAT  F.lATIcN 

Pierre  Jamet* 

D  INTRODUCTION 

Let  x  «  1R  be  the  space  variable  and  t  the  time.  Let  b^U)  and  s;(t)  be  two 
given  continuous  functions  defined  for  t  ^  0,  such  that  s^(t)  s  it)  lor  all  t. 
Let  f(x,t)  be  a  given  function  defined  for  s^ (t)  •  x  s^(t)#  t  -0  and  let  u  (x) 
be  a  given  function  defined  for  (0)  <  x  <  s^(0). 

We  want  to  solve  the  following  problem. 

Problem  iP^:  Find  a  function  u  ■  u(x,t)  such  that 

3  u  3^u 

1.1)  Lu  •  —  -  — r  *  f,  for  s  (t)  <  x  <  s  (t),  t  *  0  , 

Jt  „  .  1  2 

3x 

0 

1.2)  u(x,0)  »  u  (x),  for  s ^  ( 0 )  <  x  <  s.,10)  , 

1.3)  ulSjlt))  »  uls^lt))  »  0,  for  t  >  0  . 

We  will  assume  that  the  functions  s^(t)  and  s^tt)  satisfy  the  two  following 
hypotheses. 

1.4)  s  (t>  -  s1(t)  >  lQ  >0,  for  all  t  >  0  , 

l.M  Is^t’l  •>  s  (t)  |  ^  cQ|t'  “  *!<  f°r  aH  t  and  t*  (Lipshitz  continuity), 
with  j  *  1  or  2,  where  and  cQ  are  constants. 

The  numerical  scheme 

We  consider  a  mesh  arranged  in  the  following  way.  Let  i  be  a  space  index, 

0  1  1  .1  )«  an‘l  n  ^  0  be  a  time  index,  tfie  mesh-points  are  denoted  by  rn  ■  (x^.t11), 

where  t  "  0,  t"  <  t°+^,  x”  -  s  (tn),  xn  <  xn  ,,  x11  -  s.(tn)  (see  Figure  1). 

0  0  1  i  i+1  I  2 
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Figure  1 


We  want  to  compute  a  function  u^  which  approximates  the  solution  u  of  problem 
(P)  at  the  mesh-points  p".  We  will  denote  u"  «  u^lp").  Accordinq  to  the  method  of 
Bonnerot  and  Jamet  (11,  these  values  are  determined  by  solving  the  following  problem. 
Discrete  problem  (Ph) 


Find  (u  1,  0  <  i  <  I,  n  ^  0,  such  that 


1.6) 


n+1  n+1.  n+1 

txi+i  •  Vi,ui 


.  n  n  x  n  . 

(xi*i  -  xi-i'ui  * 


n+1 


n  +  1 


n+1 


n+1 


-  (t 


n+1 


n  uj+l  ~  U1  +  Ul+1  '  U1  _  U1  ~  Ui-1  .  “i  '  Ui~ 1 

*  n  n  n+1  n+1  n  n  n+1  n+1 

'v  -  x  x  -  x  x.  -  x  x  -  x.  ' 


i  +  1 


i+1 


i-1 
n  +  1 . 


i-1 


for  1 


i  «  1  -  1  and  n  »  0  , 


1.7) 


0 
u . 
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U°(x;>.  for 
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0,  for  n  '  0  , 
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«  I 


1.8) 


At  t  t‘l  .1 1  \  .  i.<n  t' . 


in  t  hi  ‘  >i  rn 


»  n  ,  n«  l  n«  1 .  , 

-t'ix  ,-x  . )  ♦  (x  -x  1) , 

t»l  i-l  »*1  i-l 


•  III  wi  1 1 1-  ( 1  .  i  I 


lWi  K  • 

w*u‘t<  1  J  discrete  opotator  which  appi ox i mates  the  Jifforenti.il  operator  1,  as 

wi’.l  be  shown  late*  (see  lemma  1 ) ,  and  f1'  is  a  weighted  av«tagc  of  f  (pn)  and 

i  i 

f.F"*1). 


In  this  paper,  we  will  assume  that  the  nodes  p”  are  equally  spaced  foi  each  n. 


n  n  ,  n 


l-l0)  *l4l  *  x.  -  h  -  ls,tt  )  -  s  <t  ))/I,  for  0  ‘  i  I  -  1  . 

As  already  noted  in  111,  in  the  particular  case  of  a  fixed  mesh  we  have  h*'**  •  hn  •  h 
and  «  0,  for  all  i;  hence,  the  relation  (1.6)  yields  the  classical  <"rank- 

Ni CO Ison  scheme. 

Although  successful  numerical  experiments,  as  well  as  practical  computat ions  toi 
ffxwing  boundary  problems  of  Stefan  type,  have  been  performed  by  moans  ot  the  foregoing 
method,  no  mathematical  results  had  been  proved  yet.  In  section  2  of  this  paper,  we 
prove  the  stability  of  the  method;  in  section  (,  we  establish  an  error  estimate*  and 
prove  the  convergence;  finally,  section  4  is  devoted  to  several  remarks  and  comments 
relating  the  present  method  with  other  methods  which  have  been  studied  to  solve  the 
sane  problem. 


.’)  STABILITY 


Wt*  will  u:.r  the  following  notation. 


V  I*  'l  r)l/1 

l-l 


for  any  function  *?  defined  at  the  nodes*  l*ri  for  1  *  1  l 

.  n  n*l  n 

k  •  t  -  t  . 


1  .  Alb.',  we  will  ik  1 1 


n.e* * r  «-ir  .  I 

Assume 


2.1)  kn  <_  \  min(hn,hn  for  all  n  >  0  , 

ehere  X  n  a  t*usiti  vj  n  tant .  Then,  there  exists  three  constants  *,>  md  Jr* 
which  depend  only  on  c , ,  <  and  \  such  that  : 

i  '  The  discrete  prohlew  (P^)  admits  a  univjuc  solution  -  i  u^' 1  provided 


for  all  n  >  0. 


u'  The  tolU»winu  estimate  holds 


2.2) 


lKll^->tn||u0|lh2  ♦  r •.YtVf  <ev+1  -  tv)||IvHh2) 

v»0 


for  all  n  »  0. 

Rvmaik:  because  of  assumption  (2.1),  the  condition  kn  «■  •.  is  automatical  ly  satisfied 

for  hn  small  enough. 

Proof :  For  simplicity,  we  will  omit  the  subscript  h  in  the  expression  of  the  discrete 

norm  !|  •  !lh  . 

First,  we  write  two  Inequalities  which  will  be  useful  latei .  Hypotheses  (1.4'  and 
(l.%)  imply 

2.3)  |h°*1  -  h"!  ^  ckt'min(hn,hn*1) ,  with  c  •  2c^/t^  . 

Hypotheses  (1.4),  (1.10)  and  (2.1)  imply 

i  h*l  n.  ,  n  .  n*l , 

2.4)  |x.  -  «  |  v  c’minth  ,h  ),  for  all  1  , 

with  c'  •  c^\. 

Note  that  a  condition  similar  to  (2.  3)  has  been  used  by  Lesaint  and  Raviart  I'M. 


-4- 


Now,  w«-  multiply  il, o)  by  lb*'  u^'  ♦  hnir'  ami  sub  to»  1  i  !  -  1 .  wv 


A  ♦  A  ♦  A  •  A  ,  w  X  t  h 
1  •  i  4 


,  _  v  .n  n  .  n*i  n  .n  um. 

A,  -  *  >  ih  u  -  h  u  )  (h  u  ♦  h  u  ) 

1  i  i  i  i 

v-  l 


I- 1 

l 

.  nM 
(x  .  - 

n  ,  n 
x  .  )  lu 

n*l 
♦  u 

1-1 

i  +  \ 

t«l  i*l 

i  *1 

1- 1 

y 

L 

t-i 

,  «♦» 
(xi-i  * 

ll  ,  n 
x  )  ( u 

i-l  i-l 

u*  l 

*  Vi 

l- l  ,u  t  -  u  u.  -  u.  %i 

a  -  -kn  j  l-lii - s-  -  -i - (hn+1un  ♦  hV*‘i  ♦ 

’  iii  i  hn  hu  1 


,  %  n*l  n*l  n*l  n*l 
1  - 1  u  t  -  u  u  -  u  %  \ 

lB  v  .iti"  1 - I - j-Jfcl  (hn>1u“  .  hV*1)  . 

M  hn+1  h"*1  1 


A  -  k",hn  ♦  h"*1)  l  fn(hntlu"  ♦  hV*l>  . 

4  i-1  1  1  * 

We  will  consider  separately  each  of  the  terms  A  ,A,,A  and  A  . 

1.3  4 

Estimation  of  A^:  We  have 

.  1-1  ,  ,  1-1  .  _  .  1-1 
,  ,.n  n*l  .  r  n«l,2  r  n,2i  ,  .t\+l,2  n,2,  v  n  n«l 

A  •  2h  h  I  >  (u  )  -  )  («.)  )  ♦  21  h  )  -(h))  ’  u.u 

1  '  C,  i  .  *■,  i  -.ii 

i-1  i-1  i-l 

he  cause  of  12.3),  we  have 

i  ,n*l,3  .n,2i  i.ntl  ,n  ,n  ,n*l. 

(  Ih  )  -  lh  )  I  -  [h  -  h  (h  *  h  1  - 

2.6) 

.n  .n  .  wn*l.  ,.n  ,n*l.  ,  iv  n»  1 

^  ck  (h  ♦  h  tminih  ,h  1  ^  2ck  h  h 

Hence,  hy  Cauchy's  inequality: 

A,  >  2hV+1(l  -  ok") (I  -  l)||u"n||J  ♦ 

2.7)  1  ~  " 

-  2hVM(l  ♦  ck”)  (I  -  Dllv^ll*  . 

Estimation  of  A,:  Shifting  the  index  1  m  the  second  summation  of  the  exp res 
sion  of  A^  and  taking  account  of  the  boundary  conditions  1 1.8),  we  get: 


,  1  V  ,  n+1  n  .  ,  n  n+1.  ,,n+l  n  .  n  n+1, 

A  «  -  -  l  (x  -  x  )  <u  +  u.  )  (h  u.  ♦  h  u.  )  + 


2  2  t1  '"i+1  “i+r'“i+l  T  “in' 

i»l 


1  r  ,  n+1  n,  .  n  n+1,, .n+1  n  ,n  n+1, 

2  (*i  *  Xi,(ui  +  ui  )(h  ui+l  +  h  ui+l> 


hn+l  n  n  n+1  1  n  .n+1  n  n+1,  1  ,.n+l  tn,  ,  n+1 

n  ♦  h  »  —  (h  ♦  h  )  (u^  +  )  —  —  (h  -  h  )  (ik 


.  n+1  n  n+1  n,  .n+1  .n 

(Xi+l  ’  Xi+l’  -  (xi  -  V  ’  h  ’  h 


Hence,  we  can  write 


1  ,.n+l  n  n  .n+1,  r  ,  n  n+1,,  n  n+1, 

A,  *  -  T  <h  -  h  )  h  ♦  h  )  u.  +  u.  )u.,,  +  u.  ,) 

2  4  .Mil  1  +  1  1  +  1 


1  ..n+1  n  Ir1  n+1  n  n  .  n+1  n+1 

+  -(h  *  h  )  (xi+1  -  xi+1)(ui+1  +  u.+1)(u,  -u 

-  i  (hn+1  -  h")  Y  (xn+1  -  xn)(un  +  un+1)(u"+^  -  un  ,) 
4  .‘-,1  ill  i+1  i+1 

i»l 


Applying  (2.6)  and  the  inequality 


l.  n+1  ,ni  I  n+1  ni 
I h  -  h  |  | >ci  -  x.  |  <  cc') 


which  follows  from  (2.3)  and  (2.4),  we  get 


2.8)  A, 


T  K-riK+1-i 


.  |  n  n+1 1  ■  n+1  ni  ,  ,  i  n  n+li  i  n+1 

+  c  |ui+i  +  ui+il|ui  -uil+c'  l  lui  +  ui  Nui+i 

i-l  i»l 


Apply  the  inequality 


5  2 

(a,  +  a.) (a,  +  a,)  <  )  (a.)  ,  for  all  real  numbers 

1  23  4  -  j 


we  get 


i  i«"  -  “?+1Hui+i ♦  O  ±2(I  -  +  K*in") 


"11  2  .  n  .n+1 .1  2. 


»  •  * '•»  r  •  ilwt  r^t  tht  two  eth«-i  in  the  t  i'.jht  hand  ••  \  do  mcntfu 

of  U.tfK  Hence, 

:.•»)  a,  .  a  *  euW'hi  -  ix'i^l!*’  .  |!^*l||i)  . 

t  .t  im.it  i.iii  or  a  ( 

n  n  ♦  l 

u,  u 

n  i  i  _ 

Ixit  \  •  t  Then, 

i  .n  n*l 


Men.  e , 
2 .  10) 


u  tv  n  •  l 
A  .  *-k  h  h 


1-1 


V  n  n.  n  n  , .  n 

( ( v  ,  -  V  )  -  ( V  -  V  )  )  V 

,  v  *  1  i  i  i-l  i 


i  1 


.  I-l 

n  n  n*l  r-  n  n,2 

k  h  h  iv  -  v  > 

“  1*1  1 

i“0 


.  ,  n  n  n+1  n*l  2 

. ,  1-1  ,»  ,  -  u  u  -  u  i 

.ii  nvn*l  r  I  i*l  i  .  v*l  v 

k  h  h  2  i -  ♦  - : - 


t-0  l  h" 


n*  l 


*  *  . .  II  r  (n)  n  <iv*l)  n*ln2 

h  -  k  h  h  11  -  1,!l\  \  “h 


whore  denotes  the  forward  difference  quotient  at  the  time  step  n«  i,e. 

.  n  n .  _  n 

1 1;  the  vector  with  components  tu  %  -  u  )  n  . 

»*i  i 

Fscir.ation  of  A  :  We  have 
4 

|a4i  v  xn(hn  *hn+1Mi  -  n.hnMi!<n  ♦  hnii^*1ti)rf"!i  . 

But  ,  (2.3)  implies 

.  n  .  n*l  ,  .  n.  n  .  n*l * 

h  ♦  h  <2  ♦  cfc  )min\h  ,h  )  . 

Hence 

|a4I  <_*nu  ♦  ckn)hV+1a  -  i) illujjll  ♦||^*l||)|lfnll  • 

Finally 

••.in  1a4!  nnu  *  cvn)hVMu  -  m||u£||2  ♦  ||u£4lllJ  ♦  -'ll fnll:>  • 

Taking  (2.7),  (2.9),  (2.10)  and  (2.11)  into  (2.*)  and  dividing  by  2hnh"  '(!  -  1) 
we  vjet : 


(n)  n 

X  lS, 


-7- 


,  n.n,  n*l  >2  1  .n  ,  <n>  n  .  (n«l)  n+1  2 

2.12)  a  -  I  k  (  ;uh  *  k  ■  x  uh  X  h  - 

a  .  aVMlu^ll2  *  An!linllJ.  -1th 

>in  -  c  ♦  c  (1  ♦  c1 )  *7  (2  ♦  ck"l  , 

2  4 

,n  1  ,,  .n, 

B  •  -  (2  ♦  c»  |  . 

Assume  that  the  t  ime  increments  k°  are  uniformly  bounded;  for  example,  assume  k 
Then,  we  can  replace  aIV  and  pn  in  12.12)  by  two  constants  a  •  ~  (4c  ♦  cc‘  ♦  l)  and 
B  »  l  ♦  c. 

Let  y  >  2a,  for  example  y  -  3a.  Then,  there  exists  a  constant  c  '  0  which 
depends  on  a  and  y  such  that 


1  -  ay 


-yy/2 


for  0  -  y  <  «. 


Let 


mint  2,  it  K  Then,  for  k  <  c,  we  have 


n*l  ,|  2  y  It”*1  -  tn)  | 

l“h  "  -  e 


-"II2  *  y'(tn+1  -  tB»ii?n||2  , 
n 


with  y *  •  8e^ 

The  estimate  (2.2)  follows  by  mathemat ical  induction.  The  constants  a,8,Y»*  and 
y*  depend  on  c  and  c*»  i.e.  on  co#^q  ftnd  X.  The  estimate  (2.2)  implies  the 
uniqueness  and  therefore  the  existence  of  the  solution  of  the  system  of  linear 

algebraic  equations  (1.6),  (1.7),  (1.8).  Therefore,  the  Theorem  2.1  is  proved. 


-8- 


:  :  ,V.  \*Vh  ••  t  *.  \ 


(v,^^hn  +  h"*M£-4)n+1/2 

3x  i 


.2,  n+1/2  2  n+1/2  „2  n+1/2 

+  a20(H)  -02^.  + 

3x  i  at  i 

,■1 3  n+l/2  ,3  n+1/2 

**.0(75)  . *03(7?.  * 

3x  1  3t  1 

+  terms  of  order  4,  with 

„  =  I  ,Kn+l  .n  n.2  n  2  ,  ,.n+1.2, 

a20  *  “  2  lh  ~  h  )  ((d  )  +  (h  )  +  (h  )  )  , 

_  1  ..n+1  n.  ,n.n 

a.,  *  -  7  (h  -  h  d.k  , 

11  2  1 

a  =  0  , 

02  ' 

-  -L  (v,n  .  kn+1wtii,ni2  n.  n+l  .  c  ,.n+l  2  n  2  n  1  ..n+l  n,  n,  n 
a,n  77  <h  +  h  )  (5  (h  )  -  8h  h  +  5(h  )  -  (d.)  )  d .  -  —  (h  -  h  )d.k  , 

30  12  1  1  2  1 


k"  f1  ,v.n  A  .  n+l,  ,.n  2  1  n+1  n  2  n  .n+1,  ,Ln+l  Ln.,n-i 

—  (h  +  h  )  (dj  +  -  (h  -  h  )  (h  +  h  )  +  (h  -  h  )k  ]  , 


a!2  '  °  ' 


1  n  .n+1 ,  „  n,  3 
a03  =  24  <h  +  h  >  (k  >  • 


Apply  the  inequalities 


1.  n+1  .ni  1  ,n„n  .n+1. 

I  h  -  h  I  <  -  ck  (h  +  h  ), 


1  ,n  1  ,  n  1  .  n  .n+1 , 

|d;  |  <_  c„k  <  —  c1  (h  +  h  )  , 


i1  -  0  -  2 


which  follow  from  (2.3),  (1.5)  and  (2.1).  Then,  each  of  the  coefficients  of  (3.2) 
satisfies 


| a  J  <  C'k"(hn  +hn+1)((h)2  +  (k)2)  , 
ap  — 

where  C’  is  a  positive  constant  which  depends  on  c,c'  and  cQ,  for  0  <_  a  <_  3, 

0  £  8  £  3,  a  ♦  8  s  2  or  3.  The  coefficients  of  the  terms  of  order  4  are  of  degree  4  in 
hn,  h  1  and  kn  and  satisfy  the  same  estimate  (for  which  (2.3)  is  not  even  needed). 
Hence,  after  dividing  (3.2)  by  kn (hn  +  hn+*),  we  get  (3.1)  and  the  lemma  is  proved. 
Remark :  A  Taylor  expansion  at  the  point  p"  =  (hnp"  +  hn+1p"+1)/(hn  +  hn+1)  is  more 
complicated  and  does  not  give  a  better  estimate  for  the  truncation  error. 
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xt\  ♦  1 
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•tll.l  <  .. 

0  0 


V\.  havo 


"  -  *  h"Mr»r?4ln  ,hn  .  h"M>  . 

1  l  X 


1\<  l  % 


l*t  b  f  donot«*  t  ho  *1« 


n»  1 


'  donvativo  of  t  along  t ho  vpftoi  (l  A  I  .  with 

i 


\  xl*  *'n''  '***  *  donoto  t  ho  oot  t  0*1*011*11  ng  }<o**oin1  oi*lo»  *1oi  1  v»it  ivo. 


A  Taylor  oxi'anaiou  at  t  ho  }H>mt  rn#*  *  glvo* 


i"  -  f’'*1  2  .  1  -JL  „  ,"•»  ■  ,  >  n-fu.M 

i  i  hn  ,  hnM  tft  H  V°  • 


whoio  V  is  locator  on  t  ho  straight  sogmont  (pnr  )  . 

i  i 


Hut,  taking  account  of  i».4),  wo  havo 


n  ♦  Vklfllt«.R 


Kf|  i  o  ♦  c0)\ 

Honoo,  after  using  (  X.  i)  ,  wo  got 

I .  f"*1  -I  .  .  i.Ji.i  .  i 

wtuoli  on. is  tlm  proof  of  t  tii'  lomri.i,  with 


i  1  ■  -4  “  *  C0,K  lfll,.fR  *  «»  ♦  VV'lflj.-.R  • 


.  •'  *  ort)  *  (l  »  o  )  amt  o  -  .V  f 

4  0  H  0  0  0 


Theor  om  ’.1 


tot  the  as?iiimj,'t  ions  bo  t  ti«'  ..imr  an  in  Theorem  Assume  that  oaoh  st  might  segment 

tiM 


(P^p!  1 )  ii  coiitii  ntd  m  K  foi  t 


1  -  1  and  t 


T.  1 4' t  ii  and  bo 

'  * lolutlowi  »'(  problnii  (P)  and  (P.)  rxptctivtly  and  aaavwi  u  i  C  (IH*  Hmn* 


r 


l.l>) 


|un  -  eYt  (h*  ♦  k‘)  Max  |u| 


for  0  t  T  , 


2'  n\*  4 


V*  !>.«■»  ■  I  a  «•»'!.  t  . » :  t  w!  i  h  !<  ^  •  i.  I  :»  '  ‘  0  '  • 

.  n  >n  ..  . n+1/2  ..n  +  1/2  .. 

Proof:  Let  e  *  u  -  a  .  Wo  have  (ls  u.  )  .  *  r  .#  (Lu)  *  f  •  Hence 

- —  h  n  n  n  i  l  l  i 

«■->"  '  Ik-)?  -  *?  -  ( ( U  u)  £  -  (l.u)"+,/2)  -  <f"  -  f"n/2)  -  • 

h  h  i  h  i  i  ii  i  i  i  J  r 

Applyinq  lemmas  3.1  and  3.2  and  the  inequality 

If  I  „  <  2  Max  |u|  ,  for  ii  -  1  and  2  , 

W'-'R-  21n,14  m'~'R 

(which  follows  from  (1.1)),  wo  get 

K^h-  ,  .  K1  -  C3lh^  +  k‘’  _£*  |ulm.».R  ' 

1  <.  i  <  I — 1  2<m^4 

with  *  C  f  2C,.  Then,  wo  apply  Theorem  2.1  to  the  function  e^  and  we  get 


HeXiVtV1  Max  \UX’ 
O^v  <^n- 1 

t  mm  which  the  estimate  (l.b)  of  the  theorem  follows  at  once. 


Remark  3.  1 

The  condition  that  all  straiqht  seqments  (P^P^  )  must  be  contained  in  R»  for 

1  <  i  <  I  —  1  and  tn>*  <  T,  is  satisfied  even  if  the  boundary  of  R  admits 

.  .  .  n  „  n  .  n+1.  . 

re-entrant  corners  provided  k  <_  (h  ♦  h  )/c^. 

Kcma rk  1.2:  Uniform  converqence  and  approximat ion  of  the  derivative. 

We  have 


1-1 


Max 

l<i< 1-1 


3/2, 


t^|  <  (  l  |e^|2)1/2  -  /I^T  l|e^llh-  0(h-'*).  (since  k- 


i-1 


which  implies  uniform  convergence  of  to  u.  Moreover 


n  n  , .  . ^n . 

Vl  ‘  ui  U(rui)  ~  u(Pi) 


♦  0(hl/2)  -  (p"l  +  0(h1/2)  , 
3x  i 


0  (111  1  . 


which  implies  the  uniform  converqence  of  S  u^ 


to 


3u 
3x  ' 


These  results  are  not  optimal:  the  numerical  experiments  of  |1)  have  shown 

quadratic  uniform  convergence  of  u  and  second  order  approximation  of  —  (Pn)  by 

n  Dx  i  ’ 


(dv:'7(i(u  <v2-jvi*“n" 

3u 


This  approximation  of  ~  was  used  to  compute  the  free  boundary  of  the  Stefan  problem 
with  second  order  accuracy. 

Remark  3. 3 

The  foregoing  results  can  easily  be  extended  to  the  partial  differential  operator 

r  3  .  ,  v  3u 

Lu  =  H  -  li  (a(x't) 

where  a(x,t)  is  a  given  function  of  x  and  t. 
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4  )  Cl  >MMKNTS 


Since  the  publication  of  111,  other  schemes  have  been  proposed  for  the  numerical 
solution  of  the  heat  equation  on  a  variable  mesh. 

The  scheme  of  Mon  112)  111)  is  very  closely  related  to  the  scheme  of  (1);  both 

schemes  relate  the  values  of  the  approximation  u^  at  the  six  mesh-points 

P,  ,,  P  *  P  P  !.  p"  1  and  Pn.|l  but,  the  coefficients  of  the  equations  are 

l-l  l  i*l  i-1  i  i*l  n 

different  (compare  equation  (29)  of  (1)  with  equation  (3.3)  of  (13)  with  9  *  1/2); 
both  schemes  reduce  to  the  Cranh-Nicolson  scheme  in  the  particular  case  of  a  rectan¬ 
gular  mesh.  The  scheme  of  Mori  is  derived  from  the  application  of  the  generalized 
Galerkin  method  using  time-dependent  basis  functions  as  in  Mignot  111),  whereas  our 
scheme  is  derived  from  trapezoidal  finite  elements  in  space  and  time.  Mori  proves 
the  stability  and  the  convergence  of  his  scheme  in  the  maximum  norm  under  a  severe 
condition  of  the  form  k  ^  C(hn)^,  where  C  is  a  certain  constant,  C  <  1;  under 
this  condition,  the  scheme  is  of  positive  type  and  the  maximum  principle  is  valid; 
the  proof  is  completed  for  the  Stefan  problem,  including  the  computation  of  the  free 
boundary,  in  the  case  of  one  phase  1131  and  in  the  case  of  two  phases  1 14).  Let  us 
remark  that  the  same  analysis  can  be  applied  to  our  scheme,  since  it  is  also  of  positive 
type  under  the  same  condition  as  above,  just  like  the  Crank-Nicolson  scheme.  However, 
this  condition  is  too  severe  and  is  not  satisfied  in  practical  computations. 

Another  closely  related  scheme  has  been  proposed  and  experimented  by  Miller, 

Morton  and  Baines  (101.  This  scheme  is  derived  from  the  Galerkin  method  with  time- 
dependent  basis  functions  l  (x>t)>  just  like  the  scheme  of  Mori;  but,  at  each  time 
step,  the  test  functions,  say  i|r  ,  are  different  from  the  basis  functions  and 

independent  of  the  time;  more  precisely,  il^lx.t)  »  yMxtt"  1 ) ,  for  tn  <  t  <  tn  1. 

An  integration  with  respect  to  both  variables  x  and  t  is  performed  in  the  strip 
tn  <  t  <  tn+*,  as  in  (1).  No  mathematical  analysis  of  the  method  has  been  given. 

A  general  class  of  numerical  methods  using  finite  elements  of  arbitrary  shape 
in  space  and  time  has  been  proposed  by  the  author  151  16).  These  methods  differ  from 
the  previous  ones  by  the  property  that  the  approximations  admit  a  discontinuity  with 
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i.  .t  ■  t  \  t  i*  »  »  t*  t  tn.  Uruondlt  lonal  stability  and  convenience  .u« 

•  t  v.vi  t.  i  |«-m*t.tl  (mi  at*  fi  piaf  ions  of  onirr  in  on«‘  or  »«v«r.il  *.pace- 

im-nsh'iv  .  Aj  \  1  i  it  i  ‘t»  of  thi?  met  hod  have  been  m.i.t.  with  curved  t  r4(wiotd«]  finite 
.  !.m*ni  t  r  th»  t«  l  .in  t .  »l>  1 « -m  in  one  space-d i mens  ion  (Honn**rot  and  Jam«t  IM)* 
ilv<w«v>  i  ,  in.  i  mum  rival  }u.idi  at  ur  ••  formula**  ha  v«*  been  used  to  compute  the  integrals 
involved  in  t  h«*  method,  th©  mathemat leal  results  of  |t*l  which  asium**  exact  integration 
fo  not  apply.  Tin*  nixneii  .i)  results  sh*>w  that  the  method  is  third  order  accurate.  The 
. am* •  method  has  also  k>een  applied  to  mult  i -phase  Stephan  problems  with  appearing  and 
disappearing  phase*.  141;  tn  this  cas»>,  curved  triangles  are  used  at  the  points  which 
or  respond  to  the  appearance  or  di sap)*ea ranee  of  one  of  the  phases. 

v  >t  her  methods  based  on  straight  trap©  sot  dal  finite  elements  in  s)>ace  and  time 
have  been  studied  by  Lesaint  and  Ravlart  (9]  and  l-esaint  I  Hi.  Hie  principle  of  these 
methods  is  to  write  the  heat  equation  as  a  system  of  two  differential  equations  of  the 
f  i  rst.  or  ter: 


4.  1) 


£&.»&»# 

’it  “  lx 


0 


In  19),  collocation  methods  are  proponed  for  the  numerical  solution  of  (4.1);  the' 
stability  and  convergence  of  these*  methods  are  proved  untier  the  too  stringent  condition 
that  the  domain  should  be  expanding,  i.e.  s^(t)  should  be  a  non- increasing  function 
and  r. ,  it. )  should  he  a  non-decreasing  function,  which  excludes  the  application  to  the 
Multi-phase  Stefan  problem. 

The  methods  studied  in  (H)  are  based  on  an  integration  in  the  whole  strip 
fn  <  t  <  tn+1  as  in  111.  |5),  161,  1101.  The  functions  u  and  v  are  approximated 
by  functions  and  v^.  It  is  possible  to  choose  ■  Su^/Ski  then,  the  corres¬ 

ponding  methods  are  analoqous  to  the  methods  of  111  or  (51  according  to  whether  the 
function  is  continuous  or  discontinuous  across  the  lines  t  »  t  .  But,  it  is 

also  possible  to  choose  v(i  different  from  Du^/Sx,  which  yields  methods  which  are 
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mote  vjt-ntM.il  than  the  previous  ones;  in  {articular,  the  approx  lPiat  ion  tn,, .  i.irni  f 

discont  inuitie;  with  respect  to  x.  Stability  and  convergence  proofs  ate  given, 
a-uming  exact  confutation  of  the  integrals.  These*  results  do  not  apply  to  the  method 
ot  111  stalled  in  the  present  paper,  since  this  method  was  obtained  by  using  t  lit 
trapezoidal  rule  for  the  computation  of  the  integrals. 

Finally,  let  us  mention  that  the  method  of  (1)  has  been  extended  to  the  two 
dimensional  Stefan  problem  (21,  to  systems  of  non-linear  hyperbolic  equations  with 
application  to  fluid  dynamics  |7|  and  to  problems  of  heat  diffusion  with  convection 
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